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Using the ray-optics approximation, we analyze the Casimir force in a two dimensional domain 
formed by two metallic blocks adjacent to parallel metallic sidewalls, which are separated from the 
blocks by a finite distance h. For /i > 0, the ray-optics approach is not exact because diffraction 
effects are neglected. Nevertheless, we show that ray optics is able to qualitatively reproduce a 
surprising effect recently identified in an exact numerical calculation: the force between the blocks 
varies non-monotonically with h. In this sense, the ray-optics approach captures an essential part 
of the physics of multi-body interactions in this system, unlike simpler pairwise-interaction approx- 
imations such as PFA. Furthermore, by comparison to the exact numerical results, we are able to 
quantify the impact of diffraction on Casimir forces in this geometry. 



I. INTRODUCTION 

Casimir forces, which arise from quantum vacuum 
fluctuations between uncharged surfaces [H |2l IS |4l [5], 
have attracted increasing interest in recent years due to 
rapidly improving experimental capabilities for nanoscale 
structures |3] [TJ |5]. At the same time, theoretical ef- 
forts to predict Casimir forces for geometries very un- 
like the standard case of parallel plates have begun to 
yield fruit, with several promising "exact" (arbitrary ac- 
curacy) numerical methods having been demonstrated 
for a few strong-curvature structures O HOj [TTJ [12]. In 
this paper, we explore the ability of a simple approxi- 
mate method, the ray-optics technique |T3l [H] , to bridge 
the gap between analytical calculations for simple geome- 
tries and brute-force numerics for complex structures. 
Unlike pairwise-interaction approximations such as the 
proximity-force approximation (PFA) [TS] , ray optics can 
capture multi-body interactions and thus has the poten- 
tial to predict phenomena that simpler techniques can- 
not. In particular, we show that the ray-optics approach 
can qualitatively predict a recently discovered [TTJ [T5] 
non-monotonic effect of sidewall separation on the force 
between two squares adjacent to parallel walls, as de- 
picted in Fig. [T] For a sidewall separation h — this is 
known as a "Casimir piston," and in that case has been 
has been solved exactly [17l [l^. While the ray optics 
approach is exact in this structure only for the piston 
case, its ability to capture the essential qualitative fea- 
tures for h > Q suggests a wider utility as a tool to rapidly 
evaluate different geometries in order to seek interesting 
force phenomena. Furthermore, by comparison to an ex- 
act brute-force numerical method |16j . we can evaluate 
the precise effect of diffraction (which is neglected by ray 
optics) on the Casimir force in this geometry. 

The ray optics approximation expresses the Casimir 
force in terms of a sum of contributions from all pos- 
sible classical ray paths (loops) with the same starting 
and ending point \V6\. While it is strictly valid only 



in the limit of low surface curvature, since it neglects 
diffraction effects, the rays include multiple-body inter- 
actions because there exist loops that bounce off multiple 
objects. In contrast, most other low-curvature approxi- 
mations, such as PFA [15] or other perturbative expan- 
sions [ini HO], are essentially pairwise-interaction laws, 
and can therefore miss interesting physics that occurs 
when multiple bodies are brought together. One exam- 
ple occurs in the structure depicted in Fig.[T] where there 
is a force between two square (s x s) metallic blocks sep- 
arated by a distance a that is affected by the presence of 
two infinite parallel metallic sidewalls, separated from the 
blocks by a distance h. For perfect metals in the /i = 
limit, this geometry was solved analytically in both two 
dimensions [TH] for Dirichlet boundary conditions and in 
three dimensions [211 [22] for electromagnetic fields. (By 
"two dimensions," we mean three-dimensional electro- 
magnetism restricted to z-invariant fields; equivalently, a 
combination of scalar waves with Dirichlet and Neumann 
boundary conditions, corresponding to the two polariza- 
tions.) For /i > 0, this geometry was recently solved by 
an exact computational method (that is, with no uncon- 
trolled approximations) based on numerical evaluation of 
the electromagnetic stress tensor [TTJ [16]. In this case, 
an unusual effect was observed: as h is increased from 
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FIG. 1: Schematic of a two-dimensional geometry: two metal 
squares s x s separated by a distance a, and separated from 
two adjacent metal sidewalls by a a distance h. 
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0, the force between the blocks varies non-monotonically 
with h. The attractive force between the squares actu- 
ally decreases with h down to some minimum and then 
increases toward the asymptotic limit of two isolated 
squares h oo. In the numerical solution, this non- 
monotonic effect arose as a competition between the TM 
polarization (electric field in the z direction, with Dirich- 
let boundary conditions) and the TE polarization (mag- 
netic field in the z direction, with Neumann boundary 
conditions), which have opposite dependence on h. As 
explained below, it is unclear how this non-monotonic 
effect could arise in PFA or similar methods — even if 
sidewall effects are included by restricting the pairwise 
force due to hne-of-sight interactions, it seems that the 
effect of the sidewall must always decrease monotonically 
with h. When we analyze this structure in the ray-optics 
approximation, however, we find that a similar competi- 
tion between the loops with an even and odd number of 
reflections again gives rise to a non-monotonic h depen- 
dence. 

Below, we first give a general outline of the ray-optics 
approach, explain why pairwise approximations such as 
PFA must fail qualitatively in this geometry, and then 
present our results for the structure of Fig. [T] in two di- 
mensions. This is followed by a detailed description of 
the ray-optics analysis for this structure, which involves 
a combination of analytical results for certain (even- 
reflection) paths and a numerical summation for other 
paths. 



II. RAY-OPTICS CASIMIR FORCES 

Following the framework of Ref.[T31 we express the two- 
dimensional Casimir energy via the ray-optics approxi- 
mation. The ray-optics approach recasts the Casimir en- 
ergy as the trace of the (scalar) electromagnetic Green's 
function G'(x,x') = [V^ — a;^e(x, w)]~^J(x — x'), which 
is in turn expressed as a sum over contributions from 
classical "optical" paths via saddle-point integration of 
the corresponding path-integral (this is also referred to 
as the "classical optical approximation"). The optical 
paths follow straight lines and are labeled by the number 
of specular reflections from the surfaces of the conduct- 
ing objects. In particular, the Casimir energy between 
flat surfaces for Dirichlet (77 = —1) or Neumann (r/ = 1) 
boundary conditions is given approximately by ^13j : 

Here, the length of a closed geometric path starting and 
ending at a point x is denoted by ^^(x). D,. is the set of 
points that contribute to a closed optical path reflecting 
r times from the conducting surfaces. The Casimir en- 
ergy above is thus the integral over the whole domain of 
such points. This problem reduces to computing a term- 
by-term contribution from each possible closed path, as 



determined by the specific geometrical features of the sys- 
tem under consideration. Because the Neumann (TE) 
and Dirichlet (TM) boundary conditions are given by 
the sum and difference of the even and odd paths (paths 
with even/odd numbers r of reflections), respectively, it is 
convenient to compute the contribution of even and odd 
paths separately. The Casimir force is then obtained by 
the derivative of the energy with respect to the object 
separation a. 

Equation[T|is exact for objects with zero curvature (flat 
surfaces). In the presence of curved surfaces (or sharp 
corners, which in general have measure-zero contribution 
to the set of ray paths), however, the energy will include 
additional diffractive effects that are not taken into ac- 
count by Eq. [T] One can include low-order corrections 
for small curvature [14 , but this is obviously not appli- 
cable to the case of sharp corners. There is one special 
exception, the h — "piston" : in this case, the sum over 
optical paths reduces to the method of images, which is 
exact for the interior of rectangular structures. These 
limitations are to be expected, however, since the optical 
theorem is a stationary-phase approximation. 



III. PAIRWISE-INTERACTION 
APPROXIMATIONS 

There are various pairwise-interaction force laws that 
have been proposed as approximate methods to com- 
pute Casimir forces in arbitrary geometries. The most 
well-known of these is the proximity-force approximation 
(PFA), which treats the force between two bodies as a 
pairwise sum of "parallel-plate" contributions [15]. PFA 
is exact for parallel plates, and may have low-order cor- 
rections for small curvature |23j . but is an uncontrolled 
approximation for strong curvature where it can some- 
times give qualitatively incorrect results [HJ |24l [251 EH 
E71 [5S]. Another pairwise interaction is the Casimir- 
Polder 1/r^ potential, valid in the limit of dilute media, 
which has recently been proposed as a simple (uncon- 
trolled) approximation for arbitrary geometries by renor- 
malizing it for the parallel-plate case [121 US]- In this 
section, we briefly argue why no such pairwise-interaction 
approximation can give rise to the non-monotonic depen- 
dence on h that we observe in the structure of Fig. [T] 

Of course, if one considers the pairwise interaction as a 
true two-body force, for each pair of bodies in isolation, 
then the sidewalls in the structure can have no effect 
whatsoever: the force from one sidewall on one square 
will be exactly vertical (and cancelled by the force from 
the other sidewall). However, a "lateral force" from the 
sidewalls, or equivalently an /i-dependent change in the 
attractive force between the two squares, can be obtained 
by restricting the pairwise interactions to "line of sight" 
forces. For example, when considering the force on one 
vertical edge of a square from one of the sidewalls, one 
would include contributions only from the portion of the 
sidewall that is "visible" from that edge (connected in a 
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straight line from a point on the edge to the point on 
the sidewall without passing through either square) |31j . 
For a fixed h, the hne-of-sight force on the left edge of 
a square will be different from the force on the right 
edge, since one edge will have a portion of the sidewall 
blocked by the other square, and hence there will be an 
/i-dependence of the horizontal attractive force. 

In particular, since the outside edges of the squares 
"see" (and are attracted to) a greater portion of the side- 
walls than the inside edges, the net force from the side- 
walls will always reduce the attractive force between the 
squares. Already, this contradicts the exact numerical 
calculations, in which both the Neumann force and the 
total force are greater at /i = than for /i — > cxd. 

Moreover, the effect of the sidewalls in a pairwise ap- 
proximation must always decrease with h, again contra- 
dicting our results and making non-monotonic effects im- 
possible. As h increases, two things happen: first, the 
inner edges of the square "see" a larger portion of the 
sidewalls, with area proportional to h; second, the dis- 
tance from the sidewalls to the squares increases propor- 
tional to h. The latter contribution must always domi- 
nate, however, because any pairwise force must decrease 
at least as fast as in two dimensions in order to re- 
produce the parallel-plate result. Therefore, the sidewall 
contribution must decrease monotonically at least as fast 
as in any pairwise-interaction approximation. 

Unlike pairwise-interaction approximations, we show 
below that the ray-optics approximation correctly repro- 
duces both qualitative behaviors: the total force is larger 
for h = than for /i — s- oo, and the total force is non- 
monotonic in h. 



IV. RESULTS 

Here, we present the results of our calculations for 
the general h > structure shown in Fig. [T] and com- 
pare with the numerical results from the stress-tensor 
method [TTJ [16]. It turns out that the ray-optics tech- 
nique indeed captures the non-monotonic dependence of 
the force with h, although of course the quantitative pre- 
dictions differ from the exact calculations. By defini- 
tion of ray optics, these quantitative corrections can be 
attributed to diffraction from the corners. Because we 
wish to emphasize the results of the ray-optics approach, 
rather than the details of the calculation of the differ- 
ent loop-lengths ir in Eq. [l] we defer those calculational 
details until Sec. FVland here discuss the results. 

Figures |2] and [3] show two different plots of the force vs. 
distance from the metal sidewalls h, computed via both 
Eq. [l] (solids) and the numerical stress-tensor method 
(dashed). All results are normalized by the PFA force 
between isolated squares (see captions), which are in- 
dependent of h. The bottom panel shows the con- 
tributions from Neumann boundaries (TE polarization) 
and Dirichlet boundaries (TM polarization), along with 
the total Neumann+Dirichlet force. Recall that, in the 
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FIG. 2: Casimir even (red) and odd (blue) forces vs. side- 
wall separation h normalized by the the PFA force -Fppa = 
— ftc^(3)s/87ra'^ (dashed black), computed using the ray-optics 
(solid) and stress-tensor (dashed) methods. Note that the 
ray-optics results become exact as /i ^ 0. 
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FIG. 3: Casimir force vs. sidewall separation h normalized 
by the the PFA force -Fpfa = —Hc(^{3)s/8'Ka^, computed us- 
ing the ray-optics (solid) and stress-tensor (dashed) methods. 
The Neumann (green), Dirichlet (orange) and total (black) 
forces are all normalized by the total Neumann-|-Dirichlet 
PFA force. 



ray optics approximation, the Neumann/Dirichlet forces 
are given in terms of the even- and odd-path contribu- 
tions by (even ± odd)/2, respectively, and thus the to- 
tal Neumann-|-Dirichlet force is equal to the contribu- 
tion of the even paths alone. Because the even/odd de- 
composition is more natural, in the ray optics approx- 
imation, than Neumann/Dirichlet, the top panel shows 
the even/odd contributions from the same calculations. 
(Although the stress-tensor calculation does not decom- 
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pose naturally into even and odd "reflection" contribu- 
tions, here we simply define the even/odd components as 
(Neumann ± Dirichlet) /2, respectively.) 

As h goes to zero, the ray-optics results become exact. 
The numerical computation of the stress-tensor force be- 
comes difficult for small h due to our implementation's 
uniform grid, but nevertheless the linear extrapolation of 
the numerical calculations to ft. = agree with ray optics 
to within a few percent. For h > 0, the total force for 
both the ray-optics and stress-tensor results displays a 
minimum in the range h 0.2-0.3. In particular, the 
extrema lie at ft. w 0.3 and ft ~ 0.25, respectively. Not 
only is this striking non-monotonic behavior captured by 
the ray-optics approximation, but the agreement in the 
location of the extrcmum is also excellent. 

Thus far. Figs. [2] and [3] reveal two significant differ- 
ences between the ray-optics and stress-tensor results. 
First, the forces when ft is not small differ quantitatively, 
by about 30% as ft ^ oo. Since the ray-optics approxi- 
mation is essentially obtained by dropping terms due to 
diffraction (from curved surfaces and corners) , we can at- 
tribute this quantitative difference to the diffractive con- 
tribution to the Casimir force from the finite size. In 
the large-ft limit, where the sidewalls become irrelevant 
and the ray-optics result approaches PFA, the differences 
compared to the exact solution are sometimes called edge 
effects [ini HH US] . Second, although the exact and ray- 
optics results match in the ft = limit as discussed above, 
the functional forms for small ft are quite different, at 
least for Dirichlet boundaries. In the ray optics expres- 
sions, the odd contributions have a logarithmic singular- 
ity at ft = 0, which lead to corresponding singularities 
in the Neumann and Dirichlet forces. However, in the 
exact stress-tensor calculation, only the Neumann force 
seems to display a sharp upturn in slope as ft = is ap- 
proached (although it is impossible to tell whether it is 
truly singular); the stress-tensor Dirichlet force seems to 
be approaching a constant slope (which is why we were 
able to linearly extrapolate it to ft = with good accu- 
racy) . 

Another qualitative difference appears if we look at 
the even and odd contributions in Fig. [2] whereas ray 
optics and the stress-tensor method give a similar non- 
monotonic shape for the even force, the odd forces are 
quite different. (The stress-tensor odd force is monotonic 
while the ray-optics odd force is not, while the latter goes 
to zero for large ft and the former does not.) Again, we 
attribute this to a greater sensitivity to diffraction effects, 
this time for the odd forces compared to the even forces. 
As will be argued in Sec. |V A[ the domain of integration 
and the length of even ray-optics paths have a weaker 
dependence on the corners of the squares than the odd 
paths, and thus should be less sensitive to corner-based 
diffractive effects. Fortunately, the total force depends 
solely on the even path contributions, which helps to ex- 
plain why ray-optics ultimately does effectively capture 
the non-monotonic behavior and the location of the ex- 
tremum. 



Having explored the ft-dependence of the force using 
both ray-optics and numerical stress-tensor methods, we 
now turn to Fig. |4]to study the behavior of the force as 
a function of the square separation a. Figure [4] shows 
the Casimir force vs. square separation a at constant 
h/s ^ 0.25, normalized by the PFA force between iso- 
lated squares (top panel) or by the exact force at ft = 
(bottom panel). Note that in both cases the normal- 
ization is a-dependent, unlike in Figs. [2] and [3] and any 
non-monotonicity in Fig.|4]is only an artifact of this nor- 
malization. The normalization by the PFA force allows 
us to gauge both the sidewall/edge effects (which dis- 
appear for a — > 0) and whether there is a difference in 
scaling from PFA's dependence. For the bottom 

panel, we normalize against the exact ft = force, which 
tells us whether the finite sidewall separation makes a 
difference for the large-a scaling. In both cases, we show 
a few points of the stress-tensor calculation (which be- 
came very expensive for small or large a), to get a sense 
of the accuracy of the ray-optics method at different a. 

We should expect that as a ^ 0, both the PFA and 
the ray-optics solution should approach the exact so- 
lution, because the sidewall contribution becomes neg- 
ligible. This agreement, as compared to the extrapo- 
lated numerical stress-tensor results, can be observed in 
Fig. |4](top). In contrast, for the large-a limit the ray- 
optics force appears to decay as 1 /a^ instead of for 
PFA, leading to the apparent linear growth in the top 
panel of Fig. 4 If we compare to the ft = depen- 
dence in Fig. 4 bottom), it appears to be asymptoting 
to a constant for large a, indicating that the power laws 
for ft = and ft > may be identical. However, even 
if wc had more data it would be difficult to distinguish 
the presence of, for example, logarithmic factors in this 
dependence. For the ft = case, we have analytical re- 
sults for even- and odd-path forces in Sec. |V C| from the 
analytical expressions, the odd-path ft = force clearly 
goes as and the even-path ft = force also turns 

out to have the same dependence [35]. The exact 
stress-tensor computation appears to be quite different 
from both the PFA and the ray-optics force as a function 
of a, but we were not able to go to large enough compu- 
tational cells to estimate the asymptotic power law. It is 
striking that h/s as small as 0.25 is already large enough 
to yield substantial diffraction effects in the force. 



V. DETAILS OF THE RAY-OPTICS 
COMPUTATION 

In this section, we describe the computation of the ray- 
optics approximation for the ft > squares-|-sidewalls 
structure, according to Eq. [T] This involves systemati- 
cally identifying all of the possible closed ray loops, in- 
tegrating them for a given number of reflections over the 
spatial domain, and then summing over the number of re- 
flections. It turns out that the contributions of any even 
reflection order can be integrated analytically as shown 
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FIG. 4: Casimir force vs. separation between squares a with 
constant sidewall separation h = 0.25, normalized by the 
JpFA = —ficQ(ii)sl%'Ka^ (top) and the corresponding h — Q 
force (bottom). The force is computed for both even (red) 
and odd (blue) contributions, separately. Inset: schematic 
of geometry consisting of two isolated squares with adjacent 
side walls. 



below, although the final summation over refiection order 
is still numerical. On the other hand, the integrals from 
the odd reflection orders become increasingly difficult as 
the reflection order is increased, and so we resorted to 
numerical integration for odd orders r > 5. 



A. Even paths 

Because we are dealing with perfect metals, and be- 
cause the geometry has reflection-symmetry about the x 
and y axes, it is helpful to represent the optical paths 
using an infinite periodic lattice, shown in Fig. [5] sim- 
ilax to the construction in Ref. HH The reason for this 
construction is that, because of the equal-angle law for 
specular refiections in geometric optics, a reflected ray is 
equivalent to a linearly extended ray in a mirror-image 
structure. This allows us to visualize and count the set of 



possible closed paths in a straightforward fashion. Specif- 
ically, a closed path which starts at a point x and ends 
on itself is fully determined by the set of lines that start 
at X and end in the corresponding set of image points 
on the extended lattice. The unit cell of this periodic 
construction is just two vertical black lines (of length s 
and separation a) that represent the parallel walls of the 
two squares. These are repeated with a horizontal pe- 
riod a and a vertical period s + 2h. Any path that passes 
through the gap between one of these lines and the hor- 
izontal sidewalls escapes from between the two squares 
and therefore is not counted among the closed loops. 




FIG. 5: Schematic of general 2d squares-|-sidewalls lattice. 
Lines extending from solid circles unto solid circles represent 
even reflection paths. Lines extending from open solid circles 
unto open circles represent odd reflection paths. (Here, h/a m 
0.2 and s/a « 1.) A possible even (red) and forbidden odd 
(blue) path is shown. 

In order to construct all of the closed loops that orig- 
inate at a given point x in the unit cell, we proceed as 
follows. First, we construct the mirror reflections of this 
point through the vertical lines (the boundaries of the 
squares) and the horizontal lines (the sidewalls), corre- 
sponding to reflections from these metallic walls. This 
gives us a set of points in the nearest-neighbor cells. 
Then, we construct the reflections of the nearest-neighbor 
points through their sidewalls, and so on, corresponding 
to reflections of higher and higher order. A closed loop 
is simply a line segment from the original x to one of the 
reflected points, as long as it does not escape through 
one of the gaps between the squares and sidewalls as 
explained above. Figure |5] was generated from a unit 
cell with h/a ~ 0.2 and s/a 1, and shows both even- 
reflection (solid red) and odd-reflection (dashed blue) 
paths, where in this case the odd path shown escapes 
and therefore would not be counted. In the figure, the 
points are labelled according to the number of refiections 
that generate them from the original point: solid circles 
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when the numbers of horizontal and vertical reflections 
are both even, open squares when the numbers are both 
odd, and open circles otherwise. It follows that lines con- 
necting solid circles to solid circles are even paths, and 
lines connecting solid to open circles are odd paths; the 
rays connecting solid circles and squares always escape 
and therefore do not contribute.. 

To compute the Casimir energy from these paths, the 
key quantity in Eq. [l] is the length £ of the path. Let us 
label each unit cell by (n, m) according to its horizontal 
(n) and vertical (m) offset from the cell (0, 0) where the 
original point x resides. For an even path, x must be 
connected to an even-indexed image (2n,2m), for which 
the length of the path is: 



(2na)2 + [2m {2h + s)Y 



and the angle of the path, determined by m/n, is: 



tan 0„ 



m 2h 



(2) 



(3) 



The only things left to figure out are the domain of in- 
tegration of X and the allowed (n, m) for non-escaping 
paths. If /i = 0, it is obvious that our expression reduces 
to the expression of Ref. '2T, since both the whole spatial 
domain within the unit cell and all (n, m) are allowed. 
However, when h > 0, each (n^m) will be non-escaping 
only for x in a subset of the unit cell. To determine these 
subsets, the domains of integration, we take advantage of 
the closed- loop nature of the paths to cast Eq. [TJin a dif- 
ferent light. For a given (n,m), instead of integrating 
over X and it is convenient to change variables to inte- 
grate over y and a coordinate ^ = x / cos dn.m measuring 
displacement along the path (along the 9n,m direction). 
It might seem that one should integrate ^ along the whole 
line from (0,0) to (n, m), but this may involve counting 
the same point x in the unit cell multiple times. Instead, 
to avoid over-counting loops that wrap around on them- 
selves, one instead integrates from (0,0) to (n,m) where 
h — nj gcd(?T.,m) and m = mj gcd(n,r7i) are reduced to 
lowest terms. One then obtains the following equation 
for the energy: 



he X - 

47r ^ 



«(0,0) 



£3 
n.ni 



(4) 



Although not so obvious from looking at Eq. |4] the in- 
tegral in the y direction simply counts the number of 
(n, m) paths that exist in the unit cell. Because the 
length and angle of such paths are independent of y, we 
can integrate over to obtain: 



he 
in 



2na 



(5) 



where we used J d^ = £n,m and cos9n,m = 2na/£n^m, 
between which the £n^m cancels. All that is left to figure 
out is the integral in the y direction. 

To carry the integral in the y direction, we must deter- 
mine the limits of integration, or equivalently, the range 
over which we can displace the path so that it does not 
escape. We go back to Fig. |5] for reference. Again, as 
outlined above, we extend a line from a solid circle in 
the (0,0) cell to another solid circle in the (2n, 2m)-th 
cell. For the path to be allowed, it must intercept all 
of the vertical black line segments that lie between the 
two points, i.e. at each horizontal reficction from the 
squares. Because each interception (horizontal reflec- 
tion) occurs at periodic intervals, these end up partition- 
ing the unit cell in the y direction into n sets of length 
{2h + s)/h. Note that that we divide by n, rather than 
n, because as explained above the topologically distinct 
paths are uniquely specifled by (n, m) reduced to low- 
est terms. From this simple argument, we obtain that 
the vertical displacement is {2h + s)/n — 2h, provided 
that {2h + s) /n — 2h > 0. To help visualize this result, 
it is best to think of the problem on a circle. That is, 
consider a circle of length 2h + s and partition it into 
n sets of length {2h + s)/n, as well as into two regions 
of length s and 2h. If the path is to exist, each of the 
(2h + s)/n points on the circle must not intercept the 
region marked as belonging to 2h (the air gaps between 
the squares). The result follows directly by considering 
the distance that one can displace the points before any 
of them intercepts the 2h region. 

Thus, the flnal expression for the even path energies is 
left as a sum over n and m: 



he 



n,m>0 



2h 



2h- 



2h 



Aha 



[{2naY + {2m{2h + s)Yf''^ 



(6) 



An extra factor of 2 was included in the numerator since 0, we recover the even energy expression of Ref. 17, given 
the contributions of the (— n, — m) paths are identical to in Sec. IV CI 
those of the (n, m) paths by symmetry. In the limit h = 

Although we are almost done with the even reflection 
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paths, we are missing a very important contribution to 
the energy: the PFA terms, i.e. the (n, 0) and (0,m) 
paths. The PFA energy between two parallel finite metal 
regions is a well-known result, and we include here only 
(n, 0) for completeness: 

hcsa 1 

n>0 ^ ' 

167ra2 ^'^ 

The final expression Eq. [6] must still be summed nu- 
merically over (n, to) up to some upper cutoff, but this 
is an simple computation and its convergence with the 
cutoff is discussed in Sec. IV Dl 



1. Three-reflection paths 

From Fig.|5]we verify that there are eight possible three 
reflection paths, (±2,±1) and (±1,±2) according to the 
notation of Sec. |V A[ shown in Fig. [6] The only two 
nonequivalent cases are (2, 1) and (1,2), which have dif- 
ferent domains of x integration. 
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B. Odd paths 

Unlike even reflections, odd reflection paths are quite 
tedious to compute because have less symmetry. For one 
thing, the length of a path depends not only on (n, m) 
(the offset of the unit cells being connected), but also 
on the position x of the starting point. More impor- 
tantly, the identification of the domain of x for non- 
escaping paths depends in a much more complicated way 
on (n, to), making it difficult to write down a single ex- 
pression that works for all {n,m). Therefore, we ana- 
lytically solved for the odd-path contributions only up to 
five-refiection paths, where each order requires a separate 
analysis, and treated higher-order paths by a purely nu- 
merical approach. Below, the analytical solution for the 
third-order paths is given, both to illustrate the types of 
computations that are involved and also to demonstrate 
the logarithmic singularity in the force as h —f 0. 

The results of Ref. 21] give an upper limit for the num- 
ber of odd paths r > 3 that exist in this geometry for 
h = {# paths = 2(''+3)/2). The same upper limit holds 
for h > 0, but in this case the number of paths is actu- 
ally reduced because some of the h = paths now escape. 
An analytical solution for any particular order must be- 
gin by drawing all paths for h = and then perturbing 
them for /i > to eliminate any impossible paths. At 
least for low-order paths, simple geometrical arguments 
can then determine the domain of integration. 

The coordinate dependence of odd paths arises from 
the fact that any path that reflects an odd number of 
times from the planar surfaces will be non-periodic: if we 
extend the path beyond its endpoint (= starting point), 
it will not repeat. This greatly complicates the analysis. 
For example. Fig. [5] shows one such escaping odd path. 
It turns out that as h grows, the domain of integration 
for odd paths shrinks and becomes harder to visualize. 
Moreover, if we vary h, we notice that paths for some 
{n,m), regardless of their origin x, always escape. The 
problem of determining the domain of integration and 
the allowed (n, to) seems rather difficult and does not 
seem to have a general closed-form solution. 



FIG. 6: Schematic of three-reflection paths. Blue/red repre- 
sent (1, 2)/(2, 1) paths, and the lengths £i shown are used the 
calculation of the energy. 

Inspection of Fig. [6] and Fig. [5] yields t he length of 

these paths £2.1 — 2-y/ -|- (2/i + s — y)^ . Moreover, 
we see that ^2,1 is invariant along the x-axis. For this 
particular path, the maximum displacement in the ver- 
tical direction T/max occurs when £1 sin 9 = £2 sin 9 = h, 
while the minimum Xmin and maximum a;inax horizontal 
displacements occur when £2 sin 9 = h and £1 sin 9 — h, 
respectively. Therefore, the range of integration is: 



< y < s 



ah 



< X < 



a{h + s — y) 



2h + s — y J 2h + s — y 

giving us the following expression for the energy: 

-2/1 log 



(8) 
(9) 




(10) 



In order to compute £1,2, we apply similar geometrical 
considerations. Once a gain, from Fig. [5] we obtain the 

length to be £1^2 — '^\J (a — + (2h + s)^. Similarly, 
we see that £1.2 is invariant along j/-axis, and approaches 
a minimum when x = a. Therefore, the range of integra- 
tion is: 

h< y <h + s (11) 
< X < a, (12) 

giving us the following expression for the energy: 
hcas 1 



£1,2 = - 



327r {2h + sy 



{2h + sy 



(13) 
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Adding Eq. 10 to Eq. 13 and multiplying by four to ac- 



count for the different ± sign possibilities in (n, m) yields 
the total three-reflection contribution to the energy. This 
contribution has two types of terms: a polynomial term 
from Eq. 13 and the first term of Eq. \T0\ a nd a logarith- 
mic term from the second term of Eq. |10[ The polyno- 
mial term remains at /i = 0. The more intriguing com- 
ponent of this result is the logarithmic term, which falls 
as 0{h\ogh) for small h, vanishing completely at h = 0. 
This (and similar terms for higher-order reflections) is 
the source of the logarithmic singularity in slope of the 
ray-optics odd-path Casimir force at ft- = observed in 
Fig.|2] 



2. Five-reflection paths 

The analytical solution of the five-reflection contribu- 
tion is rather complicated and is not reproduced here. 
However, it has a few interesting features that we summa- 
rize here. First, just as for the three-reflection paths, the 
five-reflection contribution has an 0{hlogh) term that 
contributes to the singularity we observe in the ray-optics 
odd-path Casimir force at /i = 0. Second, one also ob- 
tains O (log ft,) terms, which should seem to unphysically 
diverge as ft — > 0. It turns out, however, that for a path 
that gives a (log ft) contribution at small ft, there exists 
another path with a —(log ft) contribution, resulting in 
exact cancellation of any divergences. 



C. Casimir piston 

The ft = limit is the well-known "Casimir piston" ge- 
ometry. In this limit, all of the optical paths contribute to 
the Casimir energy, making it possible to compute Eq. [T] 
analytically. Ref. |21j performs this calculation in three- 
dimensions, and a similar unpublished result was ob- 
tained in two-dimensions fl7'. The two-dimensional ge- 
ometry was also solved analytically for Dirichlet bound- 
aries by another method [TB]. Here, we reproduce this 
calculation as a check on both the even- [Eq. U and odd- 



path [Eqs. 13 and 10 contributions to the energy, and 
also because the Neumann result is useful and previously 
unpublished. 

We begin by computing the ft = even-path contri- 
butions. The expression for the even energy Eq. [6] is 
much simpler than for ft > 0, because the O function 
disappears and all that is left is a polynomial function in 
terms of n, m. The even-path energy, not including the 
PFA contribution (terms where n = or m = 0, but not 
both), is given by: 



as [{naf + {msf^ 



n,m>0 



he 

'8^ 



asZ2(a, s; 3) 



(14) 
(15) 



we have identifled the summation above as the second 
order Epstein Zeta function Z2 : 



Z2(a,6;p)= J2 (("a)' + M) 

n,m>0 



2\-P/2 



(16) 



The same simplification occurs in the case of odd paths, 
the lengths of which can be found, again, by inspection 
of the lattice in Fig. [sj Again, as described in Sec. |V A| 
given a point on the unit cell x, one can determine all 
possible odd-reflection paths by drawing straight lines 
from the solid circles (x) to the open circles in the lattice. 

There are two types of open circles, each of which de- 
note two different types of paths: those that reflect across 
the X axis and have y-invariant length, as well as those 
that reflect across the y axis and have x-invariant length. 
Whichever coordinate is the invariant one gives a con- 
stant integral over the unit cell, so the double integral is 
reduced to a single integral. For example, consider those 
that reflect across the x axis and have y-invariant length: 
for these paths, we need to integrate over x in the unit 
cell and perform a double sum over (n, m). However, the 
double sum can be reduced to a single sum by eliminating 
the sum over n in favor of integrating x along the whole 
real line instead of just the unit cell. Thus, we are left 
with a single integral and a single summation. Similarly 
paths that reflect across the y axis reduce to a single in- 
tegral over y and a single sum over n. As a result of all 
this manipulation, the odd integral becomes: 



-'odd 



hes 

'32^ 



E 

m>0 



dx 



[(tos)2 + a;2]3/2 



hca -k ^ 

n>Q 



dy 



[(na)2-H?/2]3/2- 



(17) 



Again, we restrict the sum to n, m > 0, because horizon- 
tal and vertical paths are divergent terms that contribute 
only to the self-energy of the metallic walls |14j . Carrying 
out the above integrals yields the following expression for 
the odd-path energy: 



So 



hen 



dd 



(18) 



The odd-path contribution to the force is therefore oc 
1 /a^ for large a. 



D. Numerics 

To evaluate Eq. [T] numerically, we used an adaptive 
two-dimensional quadrature (cubature) algorithm [30^ to 
perform the x — {x,y) integration for each (n,m). (For 
each x and (n, m) it is easy to numerically check whether 
the path is allowed, and set the integrand to zero other- 
wise. Unfortunately, this makes the integrand discon- 
tinuous and greatly reduces the efficiency of high-order 
quadrature schemes; an adaptive trapezoidal rule might 
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have worked just as well. For very small h, this requires 
some care because the energy depends on a tiny remain- 
der between two diverging terms, as we saw in Sec. |VB 2| 
but for most h there was no difficulty. 




r (reflection order) 



FIG. 7: (Color:) Convergence error (f^+i —fr)/fr in the even 
(solid) and odd (dashed) path contributions vs. reflection 
order r for values of /i = (red), h = 0.01 (blue), and h — 0.1 
(black) (here, £r means the energy computed up to order r). 
Inset: Absolute error (fcxact —£r) /inexact for both even (solid 
black) and odd (dashed red) contributions for h = 0. 

We repeat this calculation for increasing reflection or- 
der r until the total energy £r converges to the desired 
accuracy. From general considerations, one expects the 
error |£ — f^-j for the energy from a finite r to decrease as 
0(1 /r). In particular, the path lengths £ increase propor- 
tional to r (the radius in the extended lattice), and the 
number of paths with a given length also increase propor- 
tional to r (the circumference in the extended lattice), so 
the for a given r goes as 0(l/r^). The error in 

the energy is the sum over all paths of order > r, and this 
therefore goes as 0(l/r). This scaling is verified in Fig. 7] 
which plots the relative error {£r+i —£r)/£r between the 
order-r and the order-(r+ 1) energy computations for the 
particular case of a = s = 1. In general, if the energy 
converged as 0(l/r") for some power n, one would ex- 
pect this difference to converge as 0(l/r"+-'^), and so we 
expect Fig. |7]to asymptotically go as 0(l/r^). This is 
precisely what is observed, for both even and odd paths, 
and for both h = and h > 0: all of the curves asymp- 
totically approach straight lines (on a log-log scale) with 
slope —2. 

An interesting though unfortunate result is that the 
odd-path energy requires larger r in order to obtain the 
same accuracy as the even-path energy. Though this may 
not be obvious from looking at the convergence error, it is 
clear from the inset of Fig. [7] where we plot the absolute 
error (5cxact~'S'r)/^cxact instead (at h = 0). The constant 
offset observed in the absolute errors imply that, given a 



desired accuracy constrain on the even and odd energy 
calculations, one would have to compute roughly twice 
the number of odd paths in order to obtain equivalent 
accuracy. 



VI. CONCLUDING REMARKS 

By comparing the ray-optics approximation with an 
exact brute-force calculation, we have been able to study 
both the successes and limitations of the ray-optics ap- 
proximation. On the positive side, the ray-optics ap- 
proximation is capable of capturing surprising behaviors 
that arise in closed geometries involving multiple bod- 
ies, qualitatively matching phenomena identified in exact 
brute- force calculations. In particular, the ray-optics ap- 
proximation captures the non-monotonic sidewall effects 
observed for metallic squares between parallel sidewalls, 
generalized from the Casimir piston geometry. This ef- 
fect is clearly a manifestation of the multi-body charac- 
ter of the interaction, since it does not arise in simple 
two-body force laws such as PFA. Ray optics appears to 
be unique among the current simple approximations for 
Casimir force in that it can capture such multi-body ef- 
fects, even though it cannot be quantitatively accurate in 
geometries with strong curvature. On the negative side, 
diffractive effects set in rather quickly when h is increased 
from zero, marking the agreement between ray-optics and 
the exact results only qualitative. 

This makes the ray optics approximation a promis- 
ing approach to quickly search for unusual Casimir phe- 
nomena in complicated geometries. However, since it is 
an uncontrolled approximation in the presence of strong 
curvature, any prediction by ray optics in such circum- 
stances must naturally be checked against more expen- 
sive exact calculations. There will undoubtedly be com- 
plex structures in the future where ray optics fails qual- 
itatively as well as quantitatively For instance, ray op- 
tics has more difficulty with open geometries — e.g., for 
two squares with only one sidewall, only PFA paths are 
present. On the other hand, the reach of the ray optics 
technique seems in some sense to be larger than that of 
simpler approximations such as PFA. 
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